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Abstract 

We introduce a framework for computer-aided derivation of multi-scale models. It relies on a com- 
bination of an asymptotic method used in the field of partial differential equations with term rewriting 
techniques coming from computer science. In our approach, a multi-scale model derivation is characterized 
by the features taken into account in the asymptotic analysis. Its formulation consists in a derivation of 
a reference model associated to an elementary nominal model, and in a set of transformations to apply to 
this proof until it takes into account the wanted features. In addition to the reference model proof, the 
framework includes first order rewriting principles designed for asymptotic model derivations, and second 
order rewriting principles dedicated to transformations of model derivations. We apply the method to 
generate a family of homogenized models for second order elliptic equations with periodic coefficients that 
could be posed in multi-dimensional domains, with possibly multi-domains and/or thin domains. 



1 Introduction 

There is a vast literature on multi-scale methods for partial differential equations both in applied mathe- 
matics and in many modeling areas. Among all developed methods, asymptotic methods occupy a special 
place because they have rigorous mathematical foundations and can lead to error estimates based on the 
small parameters involved in the approach. This is a valuable aspect from the model reliability point of 
view. They have been applied when a physical problem depends on one or more small parameters which 
can be some coefficients or can be related to the geometry. Their principle is to identify the asymptotic 

*This work is partially supported by the European Territorial Cooperation Programme INTERREG IV A France-Switzerland 
2007-2013. 
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model obtained when the parameters tend to zero. For instance, this method apphes in periodic homoge- 
nization, i.e. to systems consisting of a large number of periodic cells, the small parameter being the ratio 
of the cell size over the size of the complete system, see for instance |BLP78t ICD991 IJZK094] . Another 
well-developed case is when parts of a system are thin, e.g. thin plates as in [Cia], that is to say that some 
of their dimensions are small compared to others. A third kind of use is that of strongly heterogeneous 
systems e.g. |BB02j . i.e. when equation coefficients are much smaller in some parts of a system than in 
others. These three cases can be combined in many ways leading to a broad variety of configurations and 
models. In addition, it is possible to take into account several nested scales and the asymptotic charac- 
teristics can be different at each scale: thin structures to a scale, periodic structures to another, etc.... It 
is also possible to cover cases where the asymptotic phenomena happen only in certain regions or even 
are localized to the boundary. Moreover, different physical phenomena can be taken into account: heat 
transfer, solid deformations, fiuid fiow, fiuid-structure interaction or electromagnetics. In each model, the 
coefficients can be random or deterministic. Finally, different operating regimes can be considered as the 
static or the dynamic regimes, or the associated spectral problems. Today, there exists a vast literature 
covering an impressive variety of configurations. 

Asymptotic methods, considered as model reduction techniques, are very useful for complex system 
simulation and are of great interest in the software design community. They enjoy a number of advantages. 
The resulting models are generally much faster (often by several order of magnitude - depending on the kind 
of model simplification ~) to simulate than the original one and are fully parameterized. In addition, they 
do not require any long numerical calculation for building them, so they can be inserted into identification 
and optimization loops of a design process. Finally, they are of general use and they can be rigorously 
applied whenever a model depends on one or several small parameters and the error between their solution 
and nominal model solution can be estimated. 

Despite these advantages, we observe that the asymptotic modeling techniques have almost not been 
transferred in general industrial simulation software while numerical techniques, as for instance the Finite 
Element Method, have been perfectly integrated in many design tools. The main limitation factor for their 
dissemination is that each new problem requires new long hand-made calculations that may be based on 
a relatively large variety of techniques. In the literature, each published paper focus on a special case 
regarding geometry or physics, and no work is oriented in an effort to deal with a more general picture. 
Moreover, even if a large number of models combining various features have already been derived, the set of 
already addressed problems represents only a tiny fraction of those that could be derived from all possible 
feature combinations using existing techniques. 

Coming to this conclusion, we believe that what prevents the use of asymptotic methods by non- 
specialists can be formulated as a scientific problem that deserves to be posed. It is precisely the issue that 
we discuss in this paper. We would like to establish a mathematical framework for combining asymptotic 
methods of different nature and thus for producing a wide variety of models. This would allow the derivation 
of complex asymptotic models are made by computers. In this paper, we present first elements of a solution 
by combining some principles of asymptotic model derivations and rewriting methods issued from computer 
science. 

In computer science equational reasoning is usually described by rewrite rules, see |BN98j for a classical 
reference. A rewrite rule t ^ u states that every occurrence of an instance of t in a term can be replaced 
with the corresponding instance of u. Doing so, a proof based on a sequence of equality transformations is 
reduced to a series of rewrite rule applications. Rules can have further conditions and can be combined by 
specifying strategies which specify where and when to apply them, see for instance |Ter031 ICKOH ICFK051 
IBKKROli [CKLW03I 

The method developed in this paper is led by the idea of derivating models by generalization. For this 
purpose, it introduces a reference model with its derivation and a way to generate generalizations to cover 
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more cases. The level of detail in the representation of mathematical objects should be carefully chosen. 
On the one hand it should have enough precision to cover fairly wide range of models and on the other 
hand calculations should be reasonably sized. The way the generalizations are made is important so that 
they could be formulated in a single framework. 

In this paper, we select as reference problem that of the periodic homogenization of a scalar second order 
elliptic equation posed in a one-dimension domain and with Dirichlet boundary conditions. Its derivation 
is based on the use of the two-scale transform operator introduced in [ADHOOj . and reused in |BLM96j . 
We quote that homogenization of various problems using this transformation was performed according 
to different techniques in |Len97l ILen06l ILSOTI ICDOOl ICDG021 ICDG08| . Here, we follow that of [LSOT) . 
so a number of basic properties coming from this paper are stated and considered as the building blocks 
of the proofs. The complete derivation of the model is organized into seven lemmas and whose proof 
is performed by a sequence of applications of these properties. Their generalization to another problem 
requires generalization of certain properties, which is assumed to be made independently. It may also 
require changes in the path of the proof, and even adding new lemmas. The mathematical concepts are 
common in the field of partial differential equations: geometric domains, variables defined on these domains, 
functions of several variables, operators (e.g. derivatives, integrals, two-scale transform, etc.). Finally, the 
proofs of Lemmas are designed to be realizable by rewriting. 

Then, we presents a computational framework based on the theory of rewriting to express the above 
method. Each property is expressed as a rewrite rule that can be conditional, so that it can be applied or 
not according to a given logical formula. A step in a lemma proof is realized by a strategy that expresses 
how the rule applies. The complete proof of a lemma is then a sequence of such strategies. Ones we 
use have been developed in previous work |BGLj that is implemented in Maple , here we provide its 
formalization. To allow the successful application of rewriting strategies to an expression that contains 
associative and/or commutative operations, such as +, *, U, H, etc, we use the concept of rewriting modulo 
an equational theory |BN98l §11]. Without such concept one needs to duplicate the rewriting rules. 

In this work, rewriting operates on expressions whose level of abstraction accurately refiects the math- 
ematical framework. Concrete descriptions of geometric domains, functions or operators are not provided. 
Their description follows a grammar that has been defined in order that they carry enough information al- 
lowing for the design of the rewriting rules and the strategies. In some conditions of rewriting rules, the set 
of variables on which an expression depends is required. This is for example the case for the linearity prop- 
erty of the integral. Rather than introducing a typing system, which would be cumbersome and restrictive, 
we introduced a specific functionality in the form of a A-term (i.e. a program). The language of strategy 
allows this use. Put together all these concepts can express a lemma proof as a strategy, i.e. a first order 
strategy, and therefore provide a framework of symbolic computation. The concept of generalization of a 
proof is introduced as second order rewrite strategies, made with second order rewriting rules, operating on 
first order strategies. They can transform first order rewrite rules and strategies and, where appropriate, 
remove or add new ones. This framework has been implemented in the software Maple . We present 
its application to the complete proof of the reference problem and also to the generalizations of the first 
lemma, by applying second order strategies, to multi-dimensional geometrical domains, multi-dimensional 
thin domains and multi-domains. 

The paper is organized as follows. Section [2] is devoted to all mathematical aspects. This includes all 
definitions and properties, the lemmas and their proof. The principles of rewrite rules and strategies are 
formulated in Section [3l Section [4] is devoted to the theoretical framework that allows to derive a model 
and its generalizations. Implementation results are described in section [H 
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2 Skeleton of two-scale modeling 



We recall the framework of the two-scale convergence as presented in |LS07j , and the proof of the reference 
model whose implementation and extension under the form of algorithms of symbolic computation are 
discussed in Section [6l The presentation is divided into three subsections. The first one is devoted to basic 
definitions and properties, stated as Propositions. The latter are admitted without proof because they 
are assumed to be prerequisites, or building blocks, in the proofs. They are used as elementary steps in 
the two other sections detailing the proof of the convergence of the two-scale transform of a derivative, 
and the homogenized model derivation. The main statements of these two subsections are also stated as 
Propositions and their proofs are split into numbered blocks called lemmas. Each lemma is decomposed 
into steps refering to the definitions and propositions. All components of the reference model derivation, 
namely the definitions, the propositions, the lemmas and the proof steps are designed so that to be easily 
implemented and also to be generalized for more complex models. We quote that a number of elementary 
properties are used in the proof but are not explicitely stated nor cited. 



2.1 Notations, Definitions and Propositions 

Note that the functional framework used in this section is not as precise as it should be for a usual 
mathematical work. The reason is that the functional analysis is not covered by our symbolic computation. 
So, precise mathematical statements and justifications are not in the focus of this work. 

In the sequel, A C M" is a bounded open set, with measure |^|, having a "sufficiently" regular boundary 
dA and with unit outward normal denoted by ng^. We shall use the set L^{A) of integrable functions 
and the set LP{A), for any p > 0, of functions / such that S L^{A), with norm ||'y||LP(^) = (£4 \v\p 
dxY^P. The Sobolev space H^{A) is the set of functions / € L'^{A) whose gradient V/ G The 
set of p times differentiable functions on A is denoted by CP{A), where p can be any integer or 00. Its 
subset Cq{A) is composed of functions whose partial derivatives are vanishing on the boundary dA of A 
until the order p. For any integers p and q, C'^{A) C L^(A). When A = (0,ai) x ... x (0,a„) is a cuboid 
(or rectangular parallelepiped) we say that a function v defined in M" is j4-periodic if for any ^ E Z", 
v{y + Yll=i^i^i^i) — '^{y) where is the i^^ vector of the canonical basis of IR". The set of ^-periodic 
functions which are C°° is denoted by C^{A) and those which are in H^{A) is denoted by H^{A). The 
operator tr (we say trace) can be defined as the restriction operator from functions defined on the closure 
of A to functions defined on its boundary dA. Finally, we say that a sequence (ti'^)e>o € converges 
strongly in L'^{A) towards G -^^(^) when e tends to zero if lim£_j.o \ \u'' — u^\\l'^[a) — 0- The convergence 
is said to be weak if limg^o /^(^^ ~ u^)v dx = for all v S L'^{A). We write = + Os{e) (respectively 
0^(e)), where Os(e) (respectively Ow{£)) represents a sequence tending to zero strongly (respectively 
weakly) in L'^{A). Moreover, the simple notation 0(e) refers to a sequence of numbers which simply tends 
to zero. We do not detail the related usual computation rules. 

Proposition 1 [Interpretation of a weak equality] For u E L'^{A) and for any v € C^{A), 

if / u{x) v{x) dx = then u = 
J A 

in the sense of L'^{A) functions. 

Proposition 2 [Interpretation of a periodic boundary condition] For u € H^{A) and for any 

JdA 



4 



In the remainder of this section, only the dimension n = 1 is considered, the general definitions being used 
for the generalizations discussed in Section [6l 

Notation 3 [Physical and microscopic Domains] We consider an interval Q = [j Jlc'^ C M divided 

c=l 

into N{e) periodic cells (or intervals) f^c'^, of size e > 0, indexed by c, and with center x^. The translation 
and magnification (Oj''^ — Xc)/s is called the unit cell and is denoted by Q} . The variables in and in Q} 
are denoted by x'^ and x^ . 

The two-scale transform T is an operator mapping functions defined in the physical domain ^ to functions 
defined in the two-scale domain x ^ where for the reference model ^ = 0.. In the following, we shall 
denote by F, and the boundaries of 0, fi" and 0^. 

Definition 4 [Two-Scale Transform] The two-scale transform T is the linear operator defined by 

{Tu){xc,x^) = u{xc + £x^) (1) 
and then by extension T{u){x'^ , x^) = u{xc + £x^) for all x" G Qc'^ and each c in 1, .., N{e). 
Notation 5 [Measure of Domains] '^'^ = -py o,nd = JqJ^^^- 
The operator T enjoys the following properties. 

Proposition 6 [Product Rule] For two functions u, v defined in fi, 

T{uv) = {Tu){Tv). (2) 
Proposition 7 [Derivative Rule] If u and its derivative are defined in then 

i5(rn)^ 

\dx J e dx^ 

Proposition 8 [Integral Rule] If a function u G L^{^) then Tu G L^{Q'^ x Q^) and 

K° / udx = K^ [ (Tu) dx^dx\ (4) 

The next two properties are corollaries of the previous ones. 
Proposition 9 [Inner Product Rule] For two functions u, v ^ 

K° [ uv dx = K^ [ (Tu) (Tv) dx^dx^. (5) 
Proposition 10 [Norm Rule] For a function u G L^{^), 

II l|2 1 Urri ||2 /C\ 

Definition 11 [Two-Scale Convergence] A sequence G is said to be two-scale strongly (re- 

spect, weakly) convergent inL'^{^^x^^) to a limit {x^ , x^) ifTu'^ is strongly (respect, weakly) convergent 
towards u° in L'^{0 xVL^). 
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Definition 12 [Adjoint or Dual of TJ As T is a linear operator from L^(il) to L^(0^ x Q^), its adjoint 
T* is a linear operator from L'^{^1^ x 0^) to L'^(ft) defined by 

K° [ T*v u dx = [ V Tu dx^dx^. (7) 

The expression of T* can be explicited, it maps regular functions in x il^ to piecewise-constant functions 
in il. The next definition introduce an operator used as a smooth approximation of T*. 

Definition 13 [Regularization of T*] The operator B is the linear continuous operator defined from 
L2(0tt x to L'^{n) by 

Bv = v{x,^). (8) 

The nuUity condition of a function v{x'^,x^) on the boundary dQ^ x fl^ is transferred to the range Bv as 
fohows. 

Proposition 14 [Boundary Conditions of Bv] If v e C^{n^;C°°{n^)) then Bv G C^(0). 
Proposition 15 [Derivation Rule for B] If v and its partial derivatives are defined on fi" x Q} then 

dx dx^ ^dx^' ^ ' 

The next proposition states that the operator B is actuahy an approximation of the operator T* for 
0^-periodic functions. 

Proposition 16 [Approximation between T* and B] If v{x'^,x^) is continuous, continuously differ- 
entiable in and Vl^ -periodic in x^ then 

(ill 

T*v = Bv-eB{x^^)+eOs{e). (10) 

Conversely, 

Bv = T*iv) + eT*{x'p-)+eOsie). (11) 

Next, the formula of integration by parts is stated in a form compatible with the Green formula used in 
some extensions. The boundary T is composed of the two end points of the interval 0, and the unit outward 
normal nr defined on T is equal to —1 and +1 at the left- and right-endpoints respectively. 

Proposition 17 [Green Rule] Ifu,v£ H^{Q) then the traces of u and v onT are well defined and 

u— dx = [ tr{u) tr{v) nr ds{x) — [ f ^ dx. (12) 
n dx fr Jn dx 

The last proposition is stated as a building block of the homogenized model derivation. 

Proposition 18 [The linear operator associated to the Microscopic problem] For fi £ W, there 
exist 6^ G H^(^l^) solutions to the linear weak formulation 

L ""'^^ = /„. ""I^ "« ' (13) 

and is unique. Since the mapping ji i— ?■ -^—^ from M to L'^{Q^) is linear then 



Moreover, this relation can be extended to any fi G L'^{Q^). 
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2.2 Two-Scale Approximation of a Derivative 

Here we detail the reference computation of the weak two-scale limit r] = lims^oT{^) in L'^{^}^ x fl^) 
when 



Wh^n) and 



dx 



< C, (15) 
L2(n) 



C being a constant independent of e. To simplify the proof, we further assume that there exist u 







u 



1 G L2(fitt X n^) such that 



T{u^) = u° + eu^ +eOy,{e), 
i.e. 

/" {T{u') - u° - en^)?; dx^dx^ = eO(e) for ah w G L'^{n^ x rji). (16) 

We quote that Assumption (|16|) is not necessary, it is introduced to simplify the proof since it avoids some 
non-equational steps. The statement proved in the remaining of the subsection is the following. 

Proposition 19 [Two-scale Limit of a Derivative] If is a sequence bounded as in il5\) and satis- 
fying il6\) . then is independent of x^, 

=u^ -x^d^su^ (17) 

defined in x is -periodic and 
Moreover, if u'' = onT then = onVK 

The proof is split into four Lemmas corresponding to the first four blocks discussed in Section [6l the other 
three being detailed in subsection! 



Lemma 20 [First Block: Constraint on u^J is independent of x^ . 
Proof. We introduce 

n f du^ „ 

\I/ = ek}' I ——Bv dx 
Jn dx 

with V G Co°(Ott;Co°($7^)). From the Cauchy-Schwartz inequality and ([l5]), lim^^o ^' = 0. 

• Step 1. The Green formula (|12p and Proposition 1141 =^ 

^ = / dx. 

Jn dx 

• Step 2. Proposition 1151 =^ 

^ = K^ f u'B{^) dx + 0{e). 

• Step 3. Proposition [16] =^ 

^ = K^ [ u'T*{-^) dx + 0{e). 
Jn dx^ 

• Step 4. Definition [12]^ 

^ = K^[ T{u')-^ dx + 0{e). 
Jntxn^ dx' 
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Step 5. Assumption (I16p and passing to the limit when e — ?> 

k} [ dx = 0. 

Jn«xai (Jx^ 

Step 6. The Green formula ^ and = on l^tt x ^ 

1 /■ , 
K / — — j- f ax = 0. 

Step 7. Proposition [1] 

^ - 



Lemma 21 [Second Block: Two-Scale Limit of the Derivative] r] = ^p-- 
Proof. We choose v G C^{n^;C^{n^)) in 



^ = r( )7; dx^dx^. (19) 

Step 1. Definition 

Jn dx 

Step 2. Proposition [16] (to approximate T* by 5), the Green formula (jl2p . the linearity of integrals, 
and again Proposition 1161 (to approximate B by T*) 

^ = -K° / n^r*(|^)cix-^ j^u^T*{^)dx-n' j^u^T*{^^^x^)dx + 0{e). 



n 



'dx'i e 

Step 3. Definition [I2] = 



^ = -K^ [ T{u')^ dx^dx^ - - f T{u')^ dx^dx^ 



q2 



Step 4. Assumption (jl6p 



Step 5. The Green formula ()12p . Lemma [20} and passing to the limit when £ ^ 



k} j r] V dx^dx^ = I ttt^ dx^dx^ 



Step 6. Proposition [T] 
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Lemma 22 [Third Block: Microscopic Boundary Condition] is -periodic. 

Proof. In ([19]), we choose v G C^{n^;C^{n^)). 
• Step 1. The steps 1-5 of the second block =^ 



riv dx^dx^ — (u^ — x^Tr-s-)^ n-ri dx^dx^ — k,^ tt-tV dx'^dx^ = 0. 



• Step 2. Lemma [2T] =^ 

/ {u^ - x^—)v npi dxUs{x'^) = 0. (20) 

• Step 3. Definition (jl7p of and Proposition [2] 

u"*^ is 0^-periodic. (21) 

■ 

Lemma 23 [Fourth Block: Macroscopic Boundary Condition] vanishes on V'^ . 
Proof. We choose v eC^{^^), 

• Step 1. The steps 1-5 of the second block and tt^ = on T =^> 

u^'u Tiptt ds{x^)dx^ = 0. 

rttx^i 

Step 2. Proposition [1] 

u° = on r". 



2.3 Homogenized Model Derivation 

Here we provide the reference proof of the homogenized model derivation. It uses Proposition [TH] as an 
intermediary result. Let u^, the solution of a linear boundary value problem posed in 0, 

d du^ix) . 

-d^("(")^) = ^^"'' (22) 
= on r, 

where the right-hand side / € L^(il), the coefficient G C°°(0) is efi^-periodic, and there exist two 
positive constants a and /3 independent e such that 

< a < a^{x) < p. (23) 

The weak formulation is obtained by multiplication of the differential equation by a test function v G C^{il.) 
and application of the Green formula, 

K° / a'{x)^^ dx = K° [ f{x)v{x) dx. (24) 
Jn dx dx Jq 
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It is known that its unique solution is bounded as in (llSh . Moreover, we assume that for some functions 
aP{x^) and /0(x«), 

T{a') = aO and T{f) = f{x^) + 0^{e). (25) 

The next proposition states the homogenized model and is the main result of the reference proof. For 6^ a 
solution to the microscopic problem (jl3p with = 1, the homogenized coefficient and right-hand side are 
defined by 



j + dx' and f"" = J J' dx\ (26) 



Proposition 24 [Homogenized Model] The limit is solution to the weak formulation 



jjduUv^ dx«= / f^v^ dx« (27) 



Qtt dx^ dx^ 



for all gC^(17«). 

The proof is split into three lemmas. 

Lemma 25 [Fifth Block: Two-Scale Model] The couple {u^,u^) is solution to the two-scale weak 
formulation 

for any G C^(l^*) and G Cf" {VL^)) . 

Proof. We choose the test functions G C^{^^), G {VL^)) . 

• Step 1 Posing v = B{v^ + et;^) in (j24p and Proposition 1141 

St; G Co°°(f7) and / ^ed^ dB{v^ + ev') ^ f j ^(^o + ^^i) 

• Step 2 Propositions [15] and [16] 

«° i ^- (S + 1^) - »° /„ ^ ^'<"°*^^- ^ °<^>- 

• Step 3 Definition [12] and Proposition [6] 
T{a')T{^) (^ + = f T{f) dxttdx^ + 0(e). (29) 



^1 



ntixni 



Step 4 Definitions (|25p . Lemma [T9l and passing to the limit when e — )• 



which is the expected result. 



QttxQi V'^a;'' 9xV V^a;!! Qx^ J Jn»xn^ 



Lemma 26 [Sixth Block: Microscopic Problem] is solution to l[T3\) with fi = — — r and 

ax* 



dx^ dx^ dx^ 
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Proof. We choose = and v^{x^,x^) = w{x'^)ip{x^) in ^ with if G C°°{n^) and G C^i^^)- 
• Step 1 Proposition [H Lemma [20l and the Unearity of the integral =^ 



Jqi dx^ dx^ dx^ Jqi dx^ 



(30) 



Step 2 Proposition [18] with u = — — r 



as announced. 



dx^ dx^ dx^ 



Lemma 27 [Seventh Block: Macroscopic Problem] is solution to [27\ ). 

dv^ 89 
dx^ dx 



Proof. We choose G Cg°(0«) and t;i = — — G C^i^K C^{^^)) m ([28]). 



Step 1 Lemma [ 



L,^. W dx^ 8x0 Uxtt + 8x^ BxO - L^a^ ^ " • ^^^^ 



Step 2 Factorizing and definitions (|26 



^^^^ J /■ ^ ^^tt 



3 Rewriting strategies 

In this section we recah the rudiments of rewriting, namely, the definitions of terms over a signature, of 
substitution and of rewriting rules. We introduce a strategy language: its syntax and semantics in terms 
of partial functions. This language will allow us to express most of the useful rewriting strategies. 

3.1 Term, substitution and rewriting rule. 

We start with an example of rewriting rule. We define a set of rewriting variables X = {x, y} and a set 
of function symbols S = {f,g,a,b,c}. A term is a combination of elements of A" U S, for instance f{x) 
or /(a). The rewriting rule f{x) g{x) applied to a term /(a) is a two-step operation. First, it consists 
in matching the left term f{x) with the input term /(a) by matching the two occurences of the function 
symbol /, and by matching the rewriting variable x with the function symbol a. Then, the result g{a) of 
the rewriting operation is obtained by replacing the rewriting variable x occuring in the right hand side 
g{x) by the subterm a that have been associated to x. In case where a substitution is possible, as in the 
application of f{b) g{x) to /(a), we say that the rewriting rule fails. 

Definition 28 Let Ti be a countable set of function symbols, each symbol f (^T, is associated with a non- 
negative integer n, its arity ar{f) i.e. the number of arguments of f . Let X be a countable set of variables 
such that T, f] X = ^. The set of terms, denoted by T(S, X), is inductively defined by 
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• X CI T{T,,X) (i.e. every rewriting variable is a term), 

• for all f £ T, of arity n, and all ti, . . . ,t„ G T(S, A"), the expression f{ti, . . . ,tn) G T(S, Af) (i.e. the 
application of function symbols to terms gives rise to terms). 

We denote by the subset of S of the function symbols of arity n. For instance in the example / and 

g belong to Si while a and b belong to Sq. Two other common examples of terms are the expressions 

df(x) 

Integral{^l, f{x),x) and diff{f{x),x) which represent the expressions J^f{x)dx and — — . Notice that 
Integral € S3, diff^ S2, / € Si and x, € Sq. For the sake of simplicity we often keep the symbolic 
mathematical notation to express the rewriting rules. In the following we see a term as an oriented, ranked 
and rooted tree as it is usual in symbolic computation. We recall that in a ranked tree the child order is 
important. For instance the tree associated to the term Integral{Q, f{x),x) has Integral as its root which 
has three children in the order il., f, x and / has one child x. 

Definition 29 A substitution is a function a : X ^ T(S,^) such that cr{x) ^ x for x £ X. The set of 
variables that a does not map to themselves is called the domain of a, i.e. Dom[a) = {x Ci X \ cr(x) 7^ x}. If 
Dom{a) = {xi, • • • , x„} then we might write a as a = {xi ti, . . . ,Xn in} for some terms ti, ..,tn- Any 
substitution a can be extended to a mapping T(S,<Y) — >■ T{T,,X) as follows: for x X, (t{x) = cr{x), and 
for any non-variable term s = f{si, • • • , Sn), we define a{s) = f{a{si), • • • , (t(s„)). To simplify the notation 
we do not distinguish between a substitution a : X ^ 7~(S, X) and its extension a : T(S, X) T(S, X). 

The application of a substitution a to a term t, denoted by cr{t), simultaneously replaces all occurrences 
of variables in t by their a-images. 

For instance, the maping a defined by cr{x) = a is a substitution and its extension a maps f{x) and g{x) 
into f{a) and g{a). 

A rewriting rule, is a pair {l,r) where / and r are terms in 'T{T,,X); it will also be denoted by / r. 
We observe that for any two terms s,t, there exists at most one substitution a such that a{s) = t. We 
mention that a rewriting rule stands for the rule application at the top position. It is more useful to be 
able to apply a rule at arbitrary position, and more generally to specify the way rules are applied. For this 
purpose we next present a strategy language that allows to built strategies out of basic constructors. To 
this end, we introduce strategy constructor symbols ; , ®, /x, etc that do not belong to T,UX. Informally, 
the constructor " ; " stands for the composition, " © " for the left choice. Some for the application of a 
strategy to the immediate subterms of the input term, r]{x) for the fail as identity constructor, Child{j, s) 
applies the strategy s to the j^^ immediate subterm, X is a fixed-point variable, and fi is the fixed-point or 
the iterator constructor, its purpose is to define recursive strategies. For example, the strategy fxX.{s;X) 
stands for s; s; . . ., that is, it is the iteration of the application of s until a fixed-point is reached. The 
precise semantics of these constructors is given in Definition [311 

Definition 30 (Strategy) Let T be a finite set of fixed-point variables. A strategy is inductively defined 
by the following grammar: 

s ::= I r | s;s \ s0s | r]{s) \ Some{s) \ Child{j,s) \ X \ fiX.s (32) 

where j G N and X £ F. The set of strategies defined from a set of rewriting rules in T(S, X) x T(S, X) 
is denoted by S-r- 

We denote by F the failing result of a strategy and T*(S, X) = T(S, X) U F. 

Definition 31 (Semantics of a strategy) The semantics of a strategy is a function [.J : S-j-(y:,x) ~^ 
(T*(S,^) — 7> T*(S,^)) defined by its application to each grammar component: 
W(F)=F 
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F otherwise 




Isi;s2l(t) = [s2]([sil(i)) 

l[s2](t) otherwise 



(33) 



but this equation cannot be directly used to define [//X.s], since the right-hand side contains as a subphrase 
the phrase whose denotation we are trying to define. Notice that the equation ([33|) amounts to saying that 
l/iX.s] should be the least fixed-point of the operator F: 



One can show, using Knaster-Tarsky fixed-point theorem |Tar55j . that -Foo is the least fixed-point of the 
operator F, that is 



Such fixed point equations arises very often in giving denotational semantics to languages with recursive 
features, for instance the semantics of the loop "while" of the programming languages |SK95l §9, §10]. 

Example 32 Out of the basic constructors of strategies given in Definition \30[ we built up some useful 
strategies. The strategy TopDown{s) applies the strategy s to an input term t in a top down way starting 
from the root, it stops when it succeeds. That is, if the strategy s succeeds on some subterm t' of t, then it 
is not applied to the proper subterms oft'. The strategy OuterMost{s) behaves exactly like TopDown{s) 
apart that if the strategy s succeeds on some subterm t' of t, then it is also applied to the proper subterms 
oft'. The strategy BottomUp{s) (resp. InnerMost{s)) behaves like BottomUp{s) (resp. InnerMost{s)) 
but in the opposite direction, i.e. it traverses a term t starting from the leafs. The strategy Normalizer[s) 
iterates the application of s until a fixed-point is reached. The formal definition of these strategies follows: 



Let D = T*(5], X) — )• 7~*(S, X) and define !^ a partial order on D as follows: 

w Q w' iS graph{w) C graph{w'). 
Let _L be the function of empty graph, and let 

Fo = ± 

F„ = F(F„„i). 




TopDown{s) 
OuterMost{s) 
BottomUp{s) 
InnerMost{s) 
Normalizer{s) 



fiX.{s © Some{X)), 
fiX.{s; Some{X)), 
l^iX.{Some{X) © s), 
fiX.{Some{X)] s), 
^iX.{s;X). 
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Example 33 Let the variable set X = {y, z, t, and the partition S = Sq U Si U S2 of the set of function 
symbols with respect to their arity with Sq = {x, x^, x^, 50, Jl, e}, Si = {n, f , n, O, i?}, S2 = {derivative}^ 
E3 = {Integral^ with obvious definitions. We present the strategy that rewrites the expression 



= [ u{x) n{x) B{v{x^ ,x^)) dx - [ u{x) ^{B{v{x''^ ,x^))) dx + 0{e), 
Jan Jn dx 



taking into account that B{v) vanishes on the boundary 50. This term is written under mathematical form 
for simplicity, but in practice it is written from the above defined symbol of functions. Remark that the 
expression B{v(x^ ,x'^)) is a function of the variable x but this does not appear explicitly in this formulation. 
Such a case cannot appear when the grammar for terms introduced in the next section is used. We need 
the two rewriting rules 

ri := w dt w dt, 

Jan Jan 

r2 := B{v{z,y)) 0, 

and the strategy TopDown already defined. Notice that the rule ri has not effect but to detect the presence 
of the integral over the boundary. Finally, the desired strategy is: 

F := TopDown{ri]TopDown{r2)), 

and the result is 

IF\{-^) = [ u{x) n{x) B{0) dx- [ u{x) -^{B{v{x^ ,x'^))) dx + 0{e). 
Jan Jn dx 

3.2 Rewriting modulo equational theories 

So far the semantics of strategies does not take into account the properties of some function symbols, e.g. 
associativity and commutativity equahties of "+". In particular the application of the rule a + b f{a, b) 
to the term (a + c) + 6 fails. More generally we next consider the rewriting modulo an equational theory, 
i.e. a theory that is axiomatized by a set of equalities. 

For the sake of illustration, we consider the commutativity and associativity theory of +, E = {x + y = 
y + X, {x + y) + z = X + {y + z)} and the rewrite rule f{x + y) ~^ f{x) + f{y) applying the linearity rule of 
a function /. Its application to the term /((a + 6) + c) modulo E yields the set of terms {/(a + b) + f{c), 
/(a) + /(6 + c), + /(a + c)}. In the following, we define part of the semantics of a strategy modulo a 
theory, we use the notation P(T(S, X)) to denote the set of subsets of T(S, X). 

Definition 34 (Semantics of a strategy modulo) Let be E be a finitary equational theory, the seman- 
tics of a strategy modulo E is a function {.j^ : ^^(e.A') ^ (^('7'*(S, A")) ^ P(r*(S,A'))) that is partly 
defined by 

lsfi{t,,...,tn})=U^Asf{t^) 

II - rf (ti) = Uj{a,{r)} if E ^ a,{l) = t, 
Isi;s2fit) = ls2filsif{t)) 

ls,(Bs2f{t) = l^''^''^^^ ^/M(^)/{F} 
1 1^2! ^(i) otherwise 

f{t} ifisnt) = {¥} 

[|s]^(t) otherwise. 



Hs)f{t) 



The semantics of Some and Child is more complex and we do not detail it here. The semantics of the 
fixed-point operator is similar to the one given in the rewriting modulo an empty theory. 
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3.3 Conditional rewriting 



Rewriting with conditional rules, also known as conditional rewriting, extends the basic rewriting with the 
notion of condition. A conditional rewrite rule is a triplet: 

(/,r, c) 

where c is a constraint expressed in some logic. The semantics of the rule application is given by 

\ J formula aj{c) can be derived from E, 

I If otherwise. 

The set of strategies defined over rewriting rules {l,r,c) E T x T x 7^ is denoted by 3-]-,%- 

3.4 Rewriting with memory 

Some definitions or computations require storing the history of the transformations of some terms. To carry 
on, we introduce a particular function symbol M G S2 of arity two to represent the memory. Intuitively 
the term M(ti,t2) represents the term ti, besides the additional information that t2 was transformed to ti 
at an early stage. From this consideration if follows that any strategy applied to M(ii,i2) should only be 
applied to ti. Formally, we define the semantics of strategy application taking into account the memory as 
a partial function: |.]^ : S']-(^^^x) ~^ — >■ T*(S,<^)) so that: 

Isjmit) = M(|s]M(ii),t2) if t = M(ti,t2), and behaves like |.], otherwise. That is, 

Is1m(F) = F 

V(r) ifa{l) = t 
F otherwise 

Isi;S2lM(i) = Is2lM(IsilM(i)) 

'I^iImW if [siImW 7^ F 
[•S2lM(i) otherwise 



U rjuit) 

[si; S2lM(i) 



etc. 



4 A Symbolic Computation Framework for Model Deriva- 
tion 

In this section we propose a framework for the two-scale model proofs. As in Example [33l the latter 
are formulated as rewriting strategies. We notice that the following framework differs from that used in 
Example [33] in that it allows for the complete representation of the data. It does not rely on external 
structures such as hash tables. To this end, we define the syntax of the mathematical expressions by means 
of a grammar 9- 

4.1 A Grammar for Mathematical Expressions 

The grammar includes four rules to built terms for mathematical functions 3", regions Oi, mathematical 
variables V, and boundary conditions S. It involves T,Reg, T,var, ^Fun, ^Oper, and Scons which are sets 
of names of regions, variables, functions, operators, and constants so subsets of Sq. Empty expressions in 
S/jeg and T,Fun are denoted by J-^r and -Ly. The set of usual algebraic operations Sop = {+, — , x, /, is a 
subset of S2. The elements of T^Type = {Unknown, Test, Known, A^xype} C So, -^Type denoting the empty 
expression, are to specify the nature of a function, namely an unknown function (as , uP, in the proof), 
a test function (as v, v^, v^) in a weak formulation or another known function (as a^, a^, or npi). The 
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boundary conditions satisfied by a function are specified by the elements of T,bc = {d,n,pd,apd,t} C Sq 
to express tliat it satisfies Diriclilet, Neuman, periodic, anti-periodic or transmission conditions. The 
grammar also involve the symbols of functions Reg, Fun, IndexedFun, IndexedReg, IndexedVar, Oper, Var, 
and BC that define regions, mathematical functions, indexed functions or regions or variables, operators, 
mathematical variables and boundary conditions. The grammar reads as 

:J::= ®(3",^) I d I V | 

Fun(/,[V,...,V],[e,...,e],ir) I 
IndexedFun(3', V) | 

Oper{A,[7,...,3-],[V,...,V],[V,...,V],[d,...,d]) \ 

±5 I M(:r,3-), 

■Jl ::= Reg(0, [d, . . . ,d],{3l, . . . , 31, 1) \ 

IndexedReg(J, V) | 

±aj j M(3?,3?), 
V ::= Var(x,3i) | IndexedVar(V, V) | M(V,V), 
e ::= BC(c,3?,3") | M(e,e), 

where the symbols d, ®, /, K, A, x and c hold for any function symbols in Sijeg, '^Cons, ^Op, ^Fun, 
TiType, ^Oper, ^Var, and T^BC- The arguments of a region term are its region name, the list of its space 
directions (e.g. [1,3] for a plane in the variables {xi,X3)), the (possibly empty) set of subregions, the 
boundary and the outward unit normal. Those of a function term are its function name, the list of the 
mathematical variables that range over its domain, its list of boundary conditions, and its nature. Those for 
an indexed region or variable or function term are its function or variable term and its index (which should 
be discrete). For an operator term these are its name, the list of its arguments, the list of mathematical 
variable terms that it depends, the list of mathematical variable terms of its co-domain (useful e.g. for T 
when the image cannot be deduced from the initial set), and a list of parameters. Finally, the arguments of 
a boundary condition term are its type, the boundary where it applies and an imposed function if there is 
one. For example, the imposed function is set to for an homogeneous Dirichlet condition and there is no 
imposed function in a periodicity condition. We shall denote by 73j(S,0), 7j(S,0), 7v(S,0), and 7e(S,0) 
the set of terms generated by the grammar starting from the non-terminal 31, 3", V, and C. The set of all 
terms generated by the grammar (i.e. starting from 31, 3", V, or C) is denoted by 7g(S,0). Finally, we also 
define the set of terms 7g(S, Af) where each non-terminal 31, 3", V, and C can be replaced by a rewriting 
variable in X. Equivalently, it can be generated by the extension of S obtained by adding " | x" with x £ X 
in the definition of each non-terminal term. Or, by adding N ::= x, with x £ X for each non-terminal N. 

Example 35 Throughout this paper, an underlined symbol represents a shortcut whose name corresponds 
to the term name. For instance, 

Q = Reg{n, [2], 0, r, n), where T = Reg{T, [], 0, X^, 
n = Fun(n, [^'] , [] , Known) , x' = Var (x, 9!) and 9! = Reg(r2, [2] , 0, F , _L j) 

represents a region-term a one- dimensional domain named 9, oriented in the direction X2, with boundary 
F and with outward unit normal n. The shortcut F is also for a region term representing the boundary 
named F. As it can be understood from this example, except names all other fields can be void terms or 
empty lists. 

Example 36 An unknown function u{x) defined on 9 satisfying homogeneous Dirichlet boundary condition 
nix) =0 onV is represented by the function-term, 

n(x) = Fun(n, [x], Cond((i,r, 0), Unknown) where x = Var(x,il). 
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4.2 Short-cut Terms 



For the sake of conciseness, we introduce shortcut terms that are constantly used in the end of the paper: 
€ 73j(S, X), X G 7v(5], <Y) defined in / G 7dj(S, used for (discrete) indices, i G 7v(5], <Y) used as an 
index defined in /, n € 73-(S, or u{x) € 73-(S, ^) to express that it depends on the variable x and the 
indexed-term of the function u indexed by i. Similar definitions can be given for the other notations used 
in the proof as O", x", Q}, x}, fV, x^, v{x^,x}) etc. The operators necessary for the proof are the integral, 
the derivative, the two-scale transform T, its adjoint T*, and B. In addition, for some extensions of the 
reference proof we shall use the discrete sum. 

Instead of writing operator-terms as defined in the grammar, we prefer to use the usual mathematical 
expressions. The table below establishes the correspondance between the two formulations. 



I- 



du 
dx 



udx = Oper(lntegral, u, [x], [],[]), 

Oper(Partial, u, [x], [x], []), 

tr{u,x){al_) = Oper (Restriction, [x], [x^], []), 

r(n,x)(x«,xi) = Oper(r,n, [x], [x^x}], [e]), 

T*{v, [xKx}]){x) = OpeT{T*,v, [x^x}], [x], [e]), 

B{v, fe«,xi])(x) = Oper{B,v, [x^x'], \x], N), 

^Ui= Oper(Sum,n^, [i], [],[]). 
i 

The multiplication and exponentiation involving two terms / and g are written fg and as usual in 
mathematics. All these conventions have been introduced for terms in T(S,0). For terms in T(S,X) as 
those encoutered in rewriting rules, the rewriting variables can replace any of the above short cut terms. 

Example 37 The rewriting rule associated to the Green rule i f reads 

/-—V dx — [ U-— dx + [ tr(u) tr(v) n dx' . 
dx J dx J 

with the short-cuts T = Reg(r, dl, 0, -Ljj, _Lj), = Reg(Q,d2,^,T,n), x = Var(x,0) and x[ = Var(x,r). 
The other symbols u, v, x, Vt, P, dl, d2, n are rewriting variables, and for instance 

^ = Oper(Partial,n,x, [], []). 

Applying this rule according to an appropriate strategy, say the top down strategy, to a term in T(S, 0) like 

f 9f(z) 

* = y -^9{z) dz, 

for a given variable term z and function terms f,g- As expected, the result is 

-J dz + J lgndz_ 

with evident notations for n and z!_. 
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4.3 A Variable Dependency Analyzer 



The variable dependency analyzer Q is related to effect systems in computer science |MM09j . It is a 
function from TjiT., 0) to the set V{Tv{T,, 0)) of the parts of 7v(S, 0). When appHed to a term t € TjiT,, 0), 
it returns the set of mathematical variables on which t depends. The analyzer is used in the condition 
part of some rewriting rules and is inductively defined by 



e((i) = for d G Scons, 

G(x) = {x} for X G rv(S,0), 

Q{®{u,v)) = Q{u) U e{v) for u,ve 7^(S,0) and ® G Sop, 
e(±y) = 0, 

e{uix^,..,x^)) = {x^,.;x:i} for n G 7^(S,0) and x^, .., x!^ G 7^(S, 0), 

©fe) = 0(m) for u G rv(S, 0) and i G 7^(S, 0), 

0([wi, • • • ,2^!!]) = 0(2^) U • • • U 9(1/!) for ui, . . . ,7/^ G 7^(S, 0). 



eitr{u,x){x^)) = eix[), 

Q{T{u,x){xKx})) = (e(n) \ Q{x)) u e{[J,x}]) if e(x) n g(u) / 0, 

e(T*(2;, = {Q{v) \ Q{[x^,x}])) U G(x) if ei[x\x}]) D e{v) / 0, 

Q{B{v, [x^,x}]){x))) = {e{v) \ e{[xKx}])) u e{x) if e{[x\x}]) n e(i;) / 0, 



We observe that these definitions are not very general, but they are sufficient for the applications of this 
paper. To complete the definition of Q, it remains to define it on memory terms, 



Example 38 For 

Jn» Jn} ox} 
the set of mathematical variables on which ^ depends is hence inductively computed as follows: 



Q{n{x)) = {x}, @{T{u{x),x){xi,x})) = {J,x}}, @{v{xKx})) = {x^x}}, 9( ^-gl- ^ ) = {J,x^}, e{T{u{x),x) 
(xKx^) ^ ) = {x^x}}, @{J^iT{u{x),x){J,x}) ^-^l-y^ dx}) = {x«}, and e(*) = 0, that is, ^ is a 



constant function. 

4.4 Formulation of the Symbolic Framework for Model Derivation 

Now we are ready to define the framework for two-scale model derivation by rewriting. To do so, the 
rewriting rules are restricted to left and right terms G 7g(S,A') x 7g(S,,^). Their conditions c are 

formulas generated by a grammar, not explicited here, combining terms in 7g(S, X) with the usual logical 
operators in A = {V, A,], g}. It also involves operations with the dependency analyzer G. The set of terms 
generated by this grammar is denoted by 7/;(S, A", S, G, A). 



The definition of Q on the operator-terms is done case by case. 




G(n) if Q(x) C e{u) 
otherwise. 



@{M{u,v)) = G(n). 
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It remains to argue that, given a strategy s in 57-p(x;.A'),7£(E,A',g,e,A)5 the set of terms 7g(S,0) is closed 
under the apphcation of s. It is sufficient to show that for each rewriting r rule in s, the application of r 
to any term t G 7g(S,0) at any position yields a term in 7g(S,0). As an example, 7g(S,0) is not closed 
under the application of the rule x Q_, where x is a variable. But it is closed under the application of 
the linearity rule J^f + gdx f dx + g dx at any position, where f,g,x,z are rewriting variables. 

The argument is, since j\, f + gdx € 73^(S,0), then f + g £ 7j(5],0), and hence f,g G 73-(S,0). Thus, 
f dx + f^gdx £ 7j(S, 0). That is, a term in 75(S, 0) is replaced by a another term in 7j(S, 0). A more 
general setting that deals with the closure of regular languages under specific rewriting strategies can be 
found in |GGJ09| . 

A model derivation is divided into several intermediary lemmas. Each of them is intended to produce a new 
property that can be expressed as one or few rewriting rules to be applied in another part of the derivation. 
Since dynamical creation of rules is not allowed, a strategy is covering one lemma only and is operating 
with a fixed set of rewriting rules. The conversion of a result of a strategy to a new set of rewriting rules 
is done by an elementary external operation that is not a limitation for generalizations of proofs. The 
following definition summarizes the framework of symbolic computation developed in this paper. 

Definition 39 The components of the quintuplet H = (S, Af, S, ©) provide a framework for symbolic 
computation to derive multi-scale models. A two-scale model derivation is expressed as a strategy vr G 
'57^(E,A'),7i(s,A',g,e,A) ■f'^''^ which the semantics |vr]^ is applicable to an initial expression ^ G T(5],0). 

In the end of this section we argue that this framework is in the same time relatively simple, it covers the 
reference model derivation and it allows for the extensions presented in the next section. 

The grammar of terms is designed to cover all mathematical expressions occuring in the proof of the 
reference model as well as of their generalizations. A term produced by the grammar includes locally 
all useful information. This avoids the use of external tables and facilitates design of rewriting rules, 
in particular to take into account the context of subterms to be transformed. It allows also for local 
definitions, for instance a same name of variable x can be used in different parts of a same term with 
different meaning, which is useful for instance in integrals. A limitation regarding generalizations presented 
in the next section, is that the grammar must cover by anticipation all needed features. This drawback 
should be fixed in another work by supporting generalization of grammars in the same time as generalization 
of proofs. 

Each step in the proof consists in replacing parts of an expression according to a known mathematical 
property. This is well done, possibly recursively, using rewriting rules together with strategies allowing 
for precise localization. Some steps need simplifications and often use the second linearity rule of a linear 
operator, A(Xu) = XAu when A is a scalar (or is independent of the variables in the initial set of ^4). 
So variable dependency of each subterm should be determined, this is precisely what 0, the variable 
dependency analyzer, is producing. The other simplifications do not require the use of 0. In addition to 
the grammar S, the analyzer Q must be upgraded in view of each new extension. 

In all symbolic computation based on the grammar S, it is implicitely assumed that the derivatives, the 
integrals and the traces (i.e. restriction of a function to the boundary) are well defined since the regularity 
of functions is not encoded. 

Due to the algebraic nature of the mathematical proofs, this framework has been formulated by considering 
these proofs as a calculus rather than formal proofs that can be formalized and checked with a proof assistant 
|BC04irWon96j . Indeed, this is far simpler and allows, from a very small set of tools, for building significant 
mathematical derivation. To cover broader proofs, the framework must be changed by extending the 
grammar and the variable dependency analyzer only. Yet, the language Tom [BBK"'"07) does not provide a 
complete environment for the implementation of our framework since it does not support the transformation 
of rewriting rules, despite it provides a rich strategy language and a module for the specification of the 
grammar. 
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5 Transformation of Strategies as Second Order Strategies 



For a given rewriting strategy representing a model proof, one would like to transform it to obtain a 
derivation of more complex models. Transforming a strategy vr G S-j(y.,x) is achieved by applying strategies 
to the strategy vr itself. For this purpose, we consider two levels of strategies: the first order ones Sf(Y,,x) 
as defined in Definition 1301 and the strategies of second order in such a way that second order strategies can 
be applied to first order ones. That is, the second order strategies are considered as terms in a set T(S, X) 
of terms where S and X remain to be defined. Given a set of strategies Sj-(y.^x) that comes with a set of 
fixed-point variables T, we pose S D S U {-~^, ; , ®, Some, Child, r], fi} U T. Let ^ be a set of second order 
rewriting variables such that X D {X UT,) = 0. Notice that first order rewriting variables and fixed-point 
variables are considered as constants in ^(S, X), i.e. function symbols in Sq. Notice also that the arity of 
the function symbols ; ,(B, Child, fj, is two, and the arity of Some and r] is one. In particular, the rule 
I r can be viewed as the term ^ (I, r) with the symbol at the root, and the strategy iiX.s viewed as 
the term fi(X,s). This allows us to define second order strategies S^f^^-^-j by the grammar 



s:s 



s^s 



r/(s) I Some{s) | Child{j,s) \ X | ^X.s 



(34) 



Again we assume that the symbols ©,... of the second order strategies do no belong to S. The 
semantics of the strategies in S^^-^^^ are similar to the semantics of first order strategies. In addition, we 
assume that second order strategies transform first order strategies, to which they are applied, into first 
order strategies. Composing several second order strategies and applying such composition to a given first 
order strategy s provide successive transformations of s. 





Si 




Ha / 










S3 


n3 \ 


S23 


/ Ha 



Figure 1: An example of the composition of transformations of strategies. 
In the following example we illustrate the extension of an elementary strategy which is a rewriting rule. 
Example 40 For the set X = {i,j,x,x^,x^,u,e} we define si, S2, S3, and S23 four rewriting rules, 

T{—,x){x\x^) ^ ' J] — - for X en and {x\x^) e x n\ 



si 



' dx ' 



dx^ 



^,du « 1, I dT{u,x){x\x^) . ^ , / a In r^tt 

82 ■■= T{-—,x)(x*,x^) ^ J] — - — - for X en and (x\x^) G ^2" x n, 

oxi e dxl 

^,du « 1 dT(u,x)(x'^,X^) . ^ , / tt In r^a r^l 

S3 := T( — ,x)(x\x^) ^ J] — - — - for x e n^ and (x\x'-) e fi" x 0:^, 

ox £ OX^ ■' J J 



rr^^^U J 1, 1 dT(u,x)(x\x^) . ^ , / tt In r^tl r^l 

S23 ■■= T(—,x)(x\x^) -w ^ /] — - forx e n^ and e n] x $7^. 

oxi e oxf J J 
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Notice that Hi (resp. Ii2) applies the rule t^^^tt — (resp. VL-^^lj, Vt^~^VL\, and n^-^Q}-) at all of the 

OZ OZi ■'7 J 



The rule si is encountered in the reference proof, S2 is a (trivial) generalization of si in the sense that it 
applies to multi- dimensional regions Q} referenced by a set of variables {xj)i, and S3 is a second (trivial) 
generalization of si on the number of sub-regions (^j-)i ^'^^ ^> ^23 is 

a generalization combining the two previous generalizations. First, we aim at transforming the strategy si 
into the strategy S2 or the strategy S3. To this end, we introduce two second order strategies with X = {v, z} 
and S D {i, j, Q,, U^, Q,^, Partial , IndexedFun, IndexedV ar, IndexedReg} , 

^-F — ,dv dv . 

Hi := OuterMost{ — '^ — ) 

II2 ■■= OuterMost{n-^nj); Outer Mo st{n^-^Q*j); Outer Mostin^-^Q]) 

dv _ dv 
— 

dz 

positions of the input first order strategy so that 

ni(si) = S2 and n2(si) = S3. 
Once III (^nd Ii2 have been defined, they can be composed to produce S23 : 

n2ni(si) = S23 or nin2(si) = S23- 

The diagram of Figure 1 illustrates the application o/fli, II2 and of their compositions. 

The next example shows how an extension can not only change rewriting rules but also to add new 
ones. 

Example 41 To operate simplifications in the reference model, we use the strategy 

dx 

si := TopDowni— 1). 

dx 

In the generalization to multi- dimensional regions, it is replaced by two strategies involving the Kronecker 
symbol 5, usually defined as 5{i,j) = 1 if i = j and 5{i,j) = otherwise, 

(dx ' \ 
0^ ^{hj), X = yj , 

53 : = TopDown {6{i,j) —^1, i = j) , 

54 : = TopDown {6{i,j) 0, i ^ j) ■ 

The second order strategy that transforms si into the strategy N ormalizer{s2 ® S3 ® S4) is 



n := TopDown{si~^S2 ® S3 © S4). 



Notice the difference with TopDown which could not apply these rules at any position. 
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6 Implementation and Experiments 



(R) 

The framework presented in Section [4.41 has been implemented in Maple . The implementation includes, 
the language Symbtrans of strategies already presented in |BGLj . The derivation of the reference model 
presented in Section [2] has been fully implemented. It starts from an input term which is the weak 
formulation (j24p of the physical problem, 

a^lll rf^ = y f_vdx, (35) 

where a = Fun(a, [0], [ ], Known), u = Fun(n, [il], [ Dirichlet ], Unknown), v = Fun(u, [iV\, [ Dirichlet ] , Test) , 
(7 = Reg(il, [l],0,r, nn), r = Reg(r, [ ],0, ±3j, -Lj), Dirichlet = BC{Dirichlet,I^,0) and where the short- 
cuts of the operators are those of Section 14.21 The information regarding the two-scale transformation is 
provided through the test functions. For instance, in the first block the proof starts with the expression 

^ = / ^B{viJ,x^){x) dx, 
J ox 

where the test function B{v{x^,x^){x) is also an input, with v = Fun(a, i^i], [Dirichlet\i],Test), x^ = 
Var(x«,il«), xi = Var(xMli), = Reg(17tt, [1], 0, H, nf,„), ll = Reg(r», [ ±j), = Reg{n\ [1], 0, 

ri,nQi), ri = Reg(r\[ ],0,±K,-L3^), and Dirichlet^ = BC{Dirichlet]^,rl,0). 

The proof is divided into five strategies corresponding to the five blocks of the proof, each ending by 
some results transformed into rewriting rules used in the following blocks. The rewriting rules used in the 
strategies are FO-rules and can be classified into the three categories. 

• Usual mathematical rules: that represents the properties of the derivation and integration operators, 
such as the linearity, the chain rule, the Green rule, etc, 

• Specialized rules: for the properties of the two-scale calculus, as those of the two-scale transform, the 
approximation of B by the adjoint T* etc, 

• Auxiliary tools: for transformations of expressions format that are not related to operator properties 
such as the rule which transforms ipi = ip2 ™to ipi — ^2 — 0- 





Usual Rules 


Specialized Rules 


Aux. Tools 


Skeleton 


53 


14 
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Table 1: The number of first order rules used in the reference model. 



The Table [U summarizes the number of first order (FO) rules, used in the reference model, by categories. 

The reference model has been extended to cover three different kinds of configurations. To proceed to an 
extension, the new model derivation is established in a form that is as close as possible of the reference proof. 
The grammar and the dependency analyzer should be completed. Then, the initial data is determined, and 
second order (SO) strategies yielding the generalized model derivation are found and optimized. As it has 
been already mentioned, 9 and Q have already been designed to cover the three extensions. 

The first generalization is to cover multi-dimensional regions, i.e. Q C M" with n > 1. When n = 2, the 
initial term is 
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j=i j=i -> —I —I -> 

where ^ = Reg(il, [1, 2], 0,r, n^), fly = Indexed(lndexed(a, j), i), i = Var(i, /), / = Reg(/, [1, 2], 0, ±3j, ±5) 
and the choice of the test function is trivially deduced. Then, the model derivation is very similar to this 
of the reference model, see |LS07j . so much so it is obtained simply by applying the SO strategy Hidefined 
in Example 1401 This extension has been tested on the four first blocks. 

The second generalization transforms the reference model into a model with several adjacent one-dimensional 
regions (or intervals) {^k)k=i,..,m so that is still an interval i.e. $7 C M. For m = 2, the initial term is the 
same as (f35|) but with = Reg(r2, [1], {fii, ^2}-, £, nn), Qi = Reg(r2i, [1], 0, Ti, n^jj), and = Reg(Q2; [1], 
0> The two-scale geometries, all variables, all kind of functions and also the operators B and T 

are defined subregion by subregion. All definitions and properties apply for each subregion, and the proof 
steps are the same after spliting the integrals over the complete region J7 into integrals over the subregions. 
The only major change is in the fourth step where the equality u\ = at the interface between Qi and 
which is encoded as transmission conditions in the boundary conditions of and 

The third extension transforms the multi-dimensional model obtained from the first generalization to a 
model related to thin cylindrical regions, in the sense that the dimension of 57 is in the order of £ in some 
directions i £ and of the order 1 in the others i £ P e.g. Q = (0, 1) x (0, e) where = {2} and = {!}. 
The boundary T is split in two parts, the lateral part Tiat and the other parts Tother where the Dirichlet 
boundary conditions are replaced by homogeneous Neuman boundary conditions i.e. = 0. In this 
special case the integrals of the initial term are over a region which size is of the order of e so it is required 
to multiply each side of the equality by the factor 1/e to work with expressions of the order of 1. Moreover, 
the macroscopic region differs from Q, it is equal to $7" = (0, 1) when the microscopic region remains 
unchanged. In general, the definition of the adjoint T* is unchanged but {Bv){x) = v{{xi)i^jt, {x — x|)/e) 
where x| is the center of the c^^ cell in fi". It follows that the approximations (IIOI [TT]) are between T* and 
£B with y.^rs xjp. instead of V^-i x}^. With these main changes in the definitions and the preliminary 

properties, the proof steps may be kept unchanged. 





Usual Rules 


Specialized Rules 


Aux. Tools 


Mult i- D imension 


6 





4 


Thin- Region 


2 








Multi-Region 


3 









Table 2: The number of first order rules used in the three extensions. 



The mathematical formulation of the second and third extensions has been derived. This allows for the 
determination of the necessary SO-strategies, but they have not been implemented nor tested. To summa- 
rize the results about the principle of extension of strategies, we show its benefit through some statistics. 
In particular the main concerned is the reusability and the extensibility of existing strategies. The Table 
2 shows an estimate of the number of new FO-rules for the three extensions in each category and for the 
first four blocks. 

The Table 3 shows the number of SO-strategies used in each extension. Finally, the Table 4 shows, the 
ratio of the modified FO-rules and the ratio of the modified FO-strategies. The reusability ratio is high 
since most of the FO-strategies defined in the skeleton model are reused. Besides very little number of 
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Usual Rules 


Specialized Rules 


Aux. Tools 


Multi-Dimension 


9 


2 


3 


Thin-Region 











Multi-Region 


1 









Table 3: The number of second order strategies used in the extension of proofs. 



Input model 


Resulting model 


% Modi. FO-rules 


% Modi. FO-strategies 


Reference 


Multi-Dim. 


16.6% 


5% 


Multi-Dim. 


Thin 








Thin 


Multi-Reg. 





2.5% 



Table 4: The ratio of modified FO-rules and FO-strategies. 



SO-strategies is used in the extensions. This systematic way of the generation of proofs is a promising path 
that will be further validated within more complex configurations for which the proofs can not obtained 
by hand. In the future, we plan to introduce dedicated tools to aid in the design of composition of several 
extensions. 
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